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We present a first precision analysis of the transverse-momentum spectrum of gaugino pairs pro- 
duced at the LHC with center-of-mass energies of 10 or 14 TeV. Our calculation is based on a 
universal resummation formalism at next-to-leading logarithmic accuracy, which is consistently 
matched to the perturbative prediction at 0{as). Numerical results are given for the "gold-plated" 
associated production of neutrahnos and charginos for a typical benchmark point in the con- 
strained Minimal Supersymmetric Standard Model. We show that the matched resummation 
results differ considerably from the Monte Carlo predictions employed traditionally in experi- 
mental analyses. We also investigate in detail the theoretical uncertainties coming from scale and 
parton-density function variations and non-perturbative effects. 
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1. Introduction 

The Minimal Supersymmetric Standard Model (MSSM) is one of the most appealing exten- 
sions of the Standard Model of particle physics [jl|, ^. In the MSSM, the fermionic partners of the 
neutral (charged) gauge and Higgs bosons are called neutralinos (charginos) and are of particular 
importance. Stabilized by /?-symmetry, the lightest neutralino is a very promising candidate for 
the dark matter observed in the universe. Moreover, the values of gaugino masses and mixings are 
key ingredients to understand both the electroweak symmetry and the supersymmetry breakings. 
Since these particles may be light enough to be produced at current hadron colliders, they have 
been the object of particular attention, and cross sections for the production of gaugino pairs have 
been extensively studied at leading order |Q, ^, |5|, ^, ^ and next-to-leading order of perturbative 
QCD \§. 

While particle pairs are produced with zero transverse momentum (pj) in the Bom approxi- 
mation, gluon bremsstrahlung induces non-zero pj at C?(o;.v) in the strong coupling constant. The 
aim of this work is to perform an accurate calculation of the -distribution of the gaugino pairs. 
Although the use of the perturbative expansion in powers of is fully justified, when the pr of 
the produced system is of the order of its invariant mass M, the convergence of the perturbative 
expansion is spoiled by powers of large logarithmic terms, InM^/ pj, in the region where pj <^ M. 
Therefore the enhanced logarithms must be resummed to all orders in a^. 



2. Transverse-momentum resummation 



The method to systematically perform all-order resummation of classes of enhanced loga- 
rithms is well-known [^, [I^ and is performed in impact parameter (b) space. Thus we work with 
the Fourier transform W of the partonic cross section defined by 

where z = M'^/s and ^/s is the partonic center-of-mass energy. The renormalisation scale {{Xr) and 
the factorisation scale (Hf) dependencies have been removed for the sake of simplicity. After a 
Mellin transform with resptect to z, the A/^-moments of the IV-function 

W{b^,M^,N)=H{M^,N)exp[G{b^,N)] (2.2) 

factorize into a Z^-independent function H and an exponential form factor exp [G] . The function H 
encodes the full process dependence, and the form factor exp [G] resums all the terms that are order 
by order logarithmically divergent. The exponent G, which controls the soft gluon emmission, is 
thus universal, i.e. it does not depend on the final state of the process [|^]. The general expressions 
for the H- and G-functions can be found in Refs. |[l^,0]. 

After the large logarithms have been resummed in Z^-space, we have to switch back to /j^ -space 
in order to achieve a phenomenological study. Special attention has to be paid to the singularities 
in the exponent G. They are related to the presence of the Landau pole in the perturbative running 
of ttv, and a prescription is needed. In our numerical study, we follow Ref. [|l^] and deform the 
integration contour in the complex Z^-plane. Note that the divergent behavior of signals the 
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Figure 1: /77-spectra of xtX2'P^^^^ ^^e LHC. The 
LO calculation (dashed) is matched to the resummed 
calculation (full) by subtracting its fixed-order ex- 
pansion (dotted). The scale (shaded band), PDF (be- 
low) and NP (insert) uncertainties are shown, as well 
as the matched result for the LHC design energy of 
14 TeV (dot-dashed). 



Figure 2: pr-spectra of j^i^^, -pairs at the LHC. The 
matched LOh-NLL (full) and the LL (dashed) results 
are compared with the predictions of the PYTHIA 
parton shower with default (bars) and tuned (crosses) 
parameters. 



onset of non-perturbative (NP) effects at large b. These unknown effects may be extracted from 
experiment and can then be included in the exponential form factor. 

Finally, in order to conserve the full information contained in the fixed-order calculation, the 
0{as) and the resummed calculations are matched by subtracting from their sum the resummed 
cross section truncated at 0{as). 

3. Numerical results 

We now present numerical results for the associated production of neutralinos and charginos 
at the LHC with a hadronic center-of-mass energy of =10 and 14 TeV. Results for the pair 
production of neutralinos and charginos at the Tevatron and at the LHC can be found in Ref. [jlll]. 
The parton densities are evaluated in the most recent parametrisation of the CTEQ collaboration 
CTEQ6 . 6M [ pj[ ] with jUf (and jUr) set to the average mass of the final state particles, and a^^ is 
evaluated at two-loop accuracy. In the following, we choose the minimal supergravity benchmark 
point SPSla' and obtain the weak-scale supersymmetric parameters through the computer 
code SuSpect2 . 3 |[T^]. 

In Fig. we show the 777- -spectra of chargino-neutralino pairs produced at the Tevatron and at 
the LHC. As expected, the 0{as) calculation (LO) diverges at low px, but becomes finite after hav- 
ing been matched to the resummed prediction at next-to-leading logarithmic (LO-i-NLL) accuracy. 
In this region the perturbative expansion (EXP) of the resummed prediction coincides with the LO 
one. For comparison, the LO-i-NLL prediction for the 14 TeV design energy of the LHC is also pre- 



3 



Transverse-momentum resummation for gaugino-pair production at the LHC 



Jonathan DEBOVE 



sented. We also study three different sources of uncertainty: the scale variations evaluated in the 
range [?w^/2, Im^], the PDF uncertainties Apdf as defined by the CTEQ collaboration [ [l3| ] and the 
choice of three parametrisations for the NP form factor evaluated through A^p = {da^^ — do) /do 



[ ]16| , |17| ]. The scale dependence of the LO+NLL prediction is clearly improved with respect to the 
LO result and the other uncertainties are all smaller than 5% for > 5 GeV. 

In Fig. ^, we compare our LO+NLL prediction with our resummed result at leading-logarith- 
mic (LL) accuracy and two different setups for the PYTHI A6 . 4 [ |T^ ] Monte Carlo (MC) generator. 
We see that the default (STD) MC simulation is clearly improved beyond the LL approximation and 
approaches the LO-i-NLL result, but peaks at slightly smaller values of pj. This behavior can be 



improved by tuning the intrinsic pj of the partons in the hadron (AW) [11], but both simulations 
underestimate the intermediate /jj -region. 

4. Conclusion 

In summary, we have calculated the -distribution of gaugino pairs produced at the LHC at 
next-to-leading logarithmic accuracy. This has the advantage of making the perturbative predictions 
finite and reducing the scale uncertainties. We have also compared our results with a MC generator 
commonly used for experimental analyses. All these features will possibly lead to improvements 
for the experimental determination of the gaugino parameters. 
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